1/ / noise in the Two-Body Random Ensemble 
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We show that the spectral fluctuations of the Two-Body Random Ensemble (TBRE) exhibit 
1// noise. This result supports a recent conjecture stating that chaotic quantum systems are 
characterized by 1/f noise in their energy level fluctuations. After suitable individual averaging, we 
also study the distribution of the exponent a in the l// a noise for the individual members of the 
ensemble. Almost all the exponents lie inside a narrow interval around a = 1 suggesting that also 
individual members exhibit 1/f noise, provided they are individually unfolded 

PACS numbers: 05.45.Mt, 05.45.Pq, 05.65.+b, 24.60.-k 



Our understanding of quantum chaos has greatly ad- 
vanced during the last two decades. The pioneering 
work of Berry and Tabor 0, and Bohigas, Giannoni and 
Schmit snowed that there exists a close relationship 
between the energy level fluctuation properties of a quan- 
tum system and the large time scale behavior of its classi- 
cal analogue. In their seminal paper, Bohigas et at con- 
jectured that the fluctuation properties of generic quan- 
tum systems which in the classical limit are fully chaotic 
coincide with those of random matrix theory (RMT). 
This conjecture is strongly supported by experimental 
data, many numerical calculations, and analytical work 
based on semiclassical arguments. Later the interest in 
these studies was renewed with the discovery that the 
spectral statistics of quantum disordered systems is also 
well described by RMT. A review of later developments 
can be found in references [jj Q . 

Thus, RMT plays a fundamental role in quantum chaos 
studies, though it was originally introduced by Wigner 
to describe the statistical properties of high-lying en- 
ergy levels of quantum systems p|. To describe spectral 
fluctuations, RMT assumes that physical Hamiltonians 
can be substituted by "reasonable" random ensembles of 
Hamiltonian matrices, and introduces convenient statis- 
tics of their level spectra. RMT should provide the en- 
semble average of these statistics. Usually, these statis- 
tics are the nearest neighbor spacing distribution [g in- 
troduced to analyze the short range correlations, and the 
Dyson A 3 [fj statistic that allows to study long range 
correlations. 

Recently a new approach, based on traditional meth- 
ods of time series, has been proposed to analyze spectral 
fluctuations Q. In this first work we showed that the 
classical random matrix ensembles (CRME) exhibit 1/f 
noise in the fluctuations of the excitation energy. We also 
presented evidences that this is actually a universal prop- 
erty of quantum chaotic systems. The purpose of this 
brief report is to study whether the 1/f noise present in 
the spectral statistics of the CRME is also present in the 
Two-Body Random Ensemble (TBRE). As we shall see 
below, this question is very pertinent if we want to apply 
this statistic to the study of the spectral fluctuations of 



real many-body systems. 

The CRME, usually called Gaussian and Circular en- 
sembles, were chosen due to their invariant properties un- 
der certain symmetry transformations 6J. For example, 
the Gaussian Orthogonal Ensemble (GOE) is invariant 
under orthogonal transformations and is applicable to 
systems invariant under time-reversal symmetry. How- 
ever, the GOE represents systems with N-body interac- 
tions while normal systems in nature are supposed to be 
very well described by effective two-body interactions in 
the mean-field basis. TBRE was introduced to tackle this 
problem, and in this sense is more appropriated to study 
atomic nuclei, quantum dots and other mesoscopic sys- 
tems. This ensemble is constructed from a GOE in the 
2-particle Hilbert space and then propagating it to the N- 
particle Hilbert space by using the direct product structure 
of this type of spaces, (for that reason this kind of en- 
sembles are also called Embedded GOE (EGOE)) H0. 
Given the single particle states |«i >, i € 1,2, ■ ■ ■ , M, 
the two-body Hamiltonian is written as 



H = 



< v k vi\H\viVj > oj, 



(1) 



Vi<Vj,Vk<Vl 



where a\. (a Vi ) creates (destroys) a fermion in the state 
\vi >. The two body matrix elements < VkVi\H\viVj > 
are properly antisymmetryzed and are taken to be inde- 
pendent Gaussian random variables with 



< v k vi\H\viVj > = 0, 



|< v k vi\H\viVj >| = a (1 + 5 c {kl)Mj) ). 



(2) 



In this equation ~ denotes ensemble average, a is a con- 
stant and 8 C is the Kronecker delta. Then, the Hamilto- 
nian matrix in the N-particle space is defined in terms of 
these two-body matrix elements via the direct product 
structure. The only non-zero N-particle matrix elements 
are of three types 



2 



< ViV 2 ■ ■ •Wjv|i2"K«2 • • -VN > = 

< v p v 2 ■ ■ -v N \H\viv 2 ■ ■ -vn > = 

< VpV q ■ ■ -V N \H\viV 2 ■ ■ -V N > = 



^2 <v i v J \H\v l v j 

Vi<Vj <vn 

vn 

^ < v p Vi\H\viVi > 

Vi=V 2 

< V p V q \H\viV 2 >, 



> 



or those obtained by permuting the single particle states. 
All other matrix elements are zero. 

Very few analytic results are known for TBRE, con- 
trary to the classical random matrix ensembles. A very 
important result is that the level density of the TBRE 
is Gaussian in the dilute limit, which corresponds to 
(N,M) — ► oo, N/M — ► H@, instead of the semi- 
circular law for the GOE To perform a numerical 
analysis of the TBRE spectral statistics an important 
difficulty must be overcome: the TBRE is not ergodic 
[ToL [TTL . In the present context ergodicity means that 
the statistical properties of individual ensemble members 
(and hence those of the physical Hamiltonian) should al- 
ways coincide with the ensemble average. In order to 
transform TBRE into an ergodic ensemble the spectrum 
of each member must be unfolded (see the unfolding de- 
scription below) individually. In this way, GOE statistics 
is recovered. For a recent review of TBRE and more gen- 
erally EGOE see reference [l3| 

In order to establish whether 1// noise is also present 
in the spectral fluctuations of TBRE, we have studied 
four ensembles with different matrix sizes. We have 
treated N = 6 "spinless" fermions in M = 11, 12, 13 
and 14 degenerated states leading to Hilbert space di- 
mensions D — 462, 924, 1716 and 3003 respectively. The 
TBRE matrices were constructed using eqs. and 
© with cr = l. There is no relevant energy scale in 
the model and the only parameter is the dimension of 
the Hilbert space. We have diagonalized 200 matrices in 
each case to obtain the ensemble average. 

For each Hamiltonian matrix the level density p{E) 
can be separated into a smooth part ~p(E), that defines 
the main trend of the level density, and a fluctuating 
part p(E). It is well known that level fluctuations am- 
plitudes are modulated by ~p(E); therefore, to compare 
the statistical properties of different systems or different 
parts of the same spectrum, the main trend defined by 
p(E) must be removed. This procedure, called unfold- 
ing, consists in mapping the level energies Ei into new 
dimensionless levels e,-, 



E; 



E(-E»), 



1,2, 



,D, 



(4) 



We have already commented the important analytical re- 
sult of French and Wong who showed that in the TBRE 
the mean level density ~p(E) goes to Gaussian form in the 
dilute limit. However, for the dimensions of the matrices 
used in this work, the corrections to the Gaussian be- 
havior are very important and different for each matrix. 
Since the use of an accurate unfolding procedure is essen- 
(3)tial to avoid misleading results for the long range spec- 
tral correlations [lif. we have selected another method. 
Recently, the problems related to the unfolding proce- 
dure in the TBRE have been discussed deeply in Refs. 
[TH Il6| . After some tests we have chosen polynomials 
up to grade 5 to fit the accumulated level density E(-E); 
higher grades produce spurious long-range correlations. 
Finally, we have thrown 5% of the eigenvalues in the two 
spectrum edges. 

In the approach of Ref. Q , the analogy of the energy 
spectrum with a time series is established in terms of the 
4i statistic. Using the unfolded energies it is defined as 



where £(£*) is an smooth approximation to the actual 
step function S(_B) that gives the true number of energy 
levels from the ground state energy Eq and up to energy 
E. This function is given by 



(*-(*» 



(6) 



where Sj are the next-neighbor level spacings,Si = e i+1 — 
Ci, with spectral average value (s) = 1. Note that 6 q 
represents the deviation of the excitation energy of the 
(<7+l)-th unfolded level from its mean value. Moreover, it 
is closely related to the level density fluctuations. Indeed, 
we can write 



5 g — H(E q+ i) - S(£ , 9+ i) - -H(Eg + l) 



(7) 



if we appropriately shift the ground state energy; thus, it 
represents the accumulated level density fluctuations at 

E = Eq + \. 

We will profit of the formal similarity of the S q function 
with a time series to analyze its properties with numeri- 
cal techniques, normally used in the domain of complex 
systems. The most simple procedure is to study the scal- 
ing properties of its power spectrum S(k). The latter is 
defined in terms of the discrete Fourier transform 



in the usual way as 



9=1 



S(k) 



exp 



2nikq 



(8) 



(9) 



where d < D is the total number of unfolded levels con- 
sidered. In the present work d ~ 0.91? We will say that 
the spectral fluctuations of a Hamiltonian ensemble ex- 
hibit 1// noise if the ensemble averaged power spectrum 
of S q follows a power law of type 



E(E) 



p(rj)dr). 



(5) 



S(k) oc 



1 



(10) 
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with a ~ 1. For a single Hamiltonian is not clear whether 
we must impose that the bare power spectrum or some 
kind of average follows the previous power law. We shall 
explore three different possibilities below. 





(N,M)=(6,11) o 




(N,M)=(6,12) • 




(N,M)=(6,13) & 




(N,M)=(6,14) » 




iit^ Least squares fit — 
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FIG. 1: Ensemble averaged power spectra of four TBRE with 
different dimensions. The best l/f a fit is also shown. The 
curves have been displaced vertically to avoid overlapping be- 
tween them. 



Least squares fit 




log k 

FIG. 2: Three examples of power spectra for an individual 
Hamiltonian selected from the (N, M) = (6, 14) TBRE. The 
curves have also been displaced vertically to avoid overlap- 
ping. Open circles represent the raw power spectrum, open 
triangles the average values of the power spectrum obtained 
by binning the log (frequency) scale, and filled circles the 
running averaged power spectrum. The solid line represents 
the best fit in each example. The power spectrum exponents 
are a = 1.10 ± 0.07, a = 1 .10 ± 0.06, and a = 1.02 ± 0.05 
respectively. 



The results obtained for the ensemble averaged power 
spectra are shown in Fig. 2] using a log-log scale. 
It clearly seen that the calculated points spread along 
straight lines. The line slopes, i.e., the power spectrum 
exponents are obtained by means of a least-squares fit 
and their values are a = 1.09 ± 0.04, a = 1.08 ± 0.01, 
a = 1.07 ±0.01 and a = 1.07 ±0.01 for (N,M) = (6,11), 
(6,12), (6,13) and (6,14), respectively. The exponents 
are very close to one, confirming that there is 1// noise 
in the TBRE. Thus, we obtain a new and powerful check 
of the conjecture that links the spectral statistics of the 
TBRE with that of the GOE. 

Moreover, in order to study to what extent the 8 q 
statistic is also meaningful for individual spectra, we have 
randomly selected a member pertaining to the TBRE 
ensemble with (N,M) — (6,14). The upper plot of 
Fig |21 shows the power spectrum of the S q function for 
this member. Although this result suggests the existence 
of a power law, the calculated points are widely spread 
around the mean behavior, and therefore other different 
curves can be used to fit the data points. Performing 
a least square fit to a straight line we obtain a value 
a = 1.10 ± 0.07, but the error seems unreliable. Fol- 
lowing B. Mandelbrot, the problem arises because of the 
double logarithmic plot: spectral components must never 
be plotted raw, only after suitable averaging 17] . One 
of the best procedures to perform this average consists 
in dividing the high frequency portion of the logarithmic 
frequency axis into equal bins and averaging the power 
spectrum components in each bin. The result of this pro- 



cedure is shown in the second plot of fig. [5] the averaged 
data points are no more widely spread, but all of them fall 
near the mean behavior. If we perform a least square fit 
to this new set of data points we obtain a = 1.10 ± 0.06; 
therefore, the fuzzy behavior is confirmed to be related 
just to the double logarithmic plot. An alternative av- 
eraging procedure is to calculate a running or spectral 
average, (S(k)). Since in this case the dimension is large 
enough, we can divide the whole spectrum in 10 differ- 
ent sets of 256 consecutive levels. Then, the S q power 
spectrum is calculated for each level set and in order to 
reduce fluctuations and clarify the main trend a running 
average is performed using these sets. The bottom plot 
of Fig. [^displays the result of this calculation using open 
circles. A least squares fit leads to the following exponent 
a = 1.02 ±0.05, which is very similar to that obtained by 
means of an ensemble average. This result is compatible 
with the ergodic properties of TBRE once the spectra are 
individually unfolded. Note that the different lengths of 
the three power spectra shown in this figure are due to 
the fact that in the first and the second cases we are ac- 
tually using the whole sequence, but in the third one we 
use sequences of 256 consecutive levels. 

To make this discussion more quantitative, we have 
calculated the exponent a for the 200 matrices of the en- 
semble with (N, M) = (6, 14) using the binning method 
previously described. The average value is < a >= 1.06 
and the width of the distribution is a a — 0.08. Figure 
13 shows an histogram of the distribution together with a 
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FIG. 3: Histogram of the distribution of a for the (N, M) = 
(6, 14) TBRE. 



Gaussian defined by the previous parameters that seems 
to ht the data very well. Although this result has been 
obtained for a particular sample of a particular ensemble 



((N,M) = (6, 14)), it suggests that individual members 
also are characterized by 1// noise. 

We have confirmed that the spectral fluctuations of 
the TBRE exhibit l/f noise. This behavior supports 
the previously stated conjecture that chaotic quantum 
systems are characterized by l/f noise in their energy 
level fluctuations. We have also shown that individual 
members have l/f noise in their excitation energy fluc- 
tuations provided they are individually unfolded and the 
power spectrum of the S q function is appropriately av- 
eraged. Actually, the distribution of the a exponent in 
the 1// Q law is a Gaussian centered near a = 1 with a 
quite small width. Therefore, the spectral fluctuations 
of atomic nuclei, quantum dots and mesoscopic systems 
can be studied by means of the scaling properties of the 
power spectrum of the S q function. The advantages of 
this new statistic are perfectly used in some recent works 
about the nuclear masses El EH, where the l/f noise 
of different series of fluctuations in the nuclear masses 
along the nuclear chart was explored. Depending of the 
definition of the fluctuations, different exponents in the 
power law were found, l/f noise was shown to be a pow- 
erful tool to investigate spectral correlations in this kind 
of experimental data. 
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